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Abstract. Resummation of the soft gluon radiative corrections for the quark-vector boson vertex is per- 
formed within the path-integral (world-line) approach. The leading order expression for the vacuum aver- 
aged Wilson integral for an arbitrary gauge field is found in n-dimensional space-time. The cusp anomalous 
dimension of the color non-singlet Sudakov form factor of on-mass-shell quark is calculated in an arbitrary 
covariant gauge in the one-loop order, and the leading double-logarithmic asymptotical behavior is obtained 
from the corresponding evolution equation. 

PACS. 12.38.Cy Summation of perturbation theory - 12.38.Lg Other nonperturbative calculations 



1 Introduction 

Resummation of the logarithmic corrections in high-energy 
processes is required to justify effectiveness of the Stan- 
dard Model at highest energies accessible at modern and 
future colliders (see, e.g., [1] and references therein). In 
electroweak processes, the Sudakov corrections can influ- 
ence profoundly the cross sections at e + e~ linear colliders 
at TeV energies, and the precise evaluation of them is 
quite important in search for the New Physics, as well as 
to test the predictive power of the Standard Model [2] . In 
the strong interactions sector, the form factors of quarks 
are the most elementary entities exhibiting the double- 
logarithmic behavior. The quark form factors being, of 
course, unobservable quantities, enter into the quark-photon 
and quark-gluon vertices in the calculations of various 
QCD processes at the partonic level, and are of a spe- 
cial theoretical as well as phenomcnological importance 
[3,4]. Very recently, the investigations of the form factors 
of quarks has been connected with the progress in experi- 
mental study of the constituent quarks and search for the 
perturbative and nonperturbative effects in their structure 
from low to high energy [5-7] . 

The form factors of the elementary fermions at large 
momenta transfers in gauge theories — in QED, and later 
in QCD — had been studied extensively since the fifties 
[8-10] [11,12]. In the color singlet case — corresponding to 
the elastic ora-shell quark scattering in an external elec- 



tromagnetic field — the exponentiation of the infrared sin- 
gularities has been proved and the correct leading asymp- 
totic at high energy has been obtained [10,15,16,18]: the 
high-energy behavior of the on-mass-shell quark form fac- 
tor can be described in terms of the perturbative evolution 
equation 

pvrr [q2] jpmi [q2] = 

expf- [ Q p ln^- r cusp [a s (fj,)]+NLo] . (1) 
V JQl 2^ A* J 

In this case, the so-called cusp anomalous dimension reads 
in one-loop order [16] 



r^Ua s )^-C F + 0{al)>Q 



(2) 
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and, therefore, the IR sensitive part of the electromagnetic 
quark form factor experiences the Sudakov suppression at 
large Q 2 . On the other hand, it has been shown recently, 
that in the similar electroweak reactions the cusp anoma- 
lous dimension has the opposite sign for charged bosons 
, and then the enhancement takes place instead of the 
suppression [19]. 

Here we apply the powerful Wilson integrals techniques 
[13, 16-18] to study of the Sudakov resummation of soft 
gluon radiative corrections to the quark-vector boson ver- 
tex with large transferred momentum. In particular, we 
concentrate on the ggg-vertex, since qq'y-vevtex. can be 
trivially restored from the latter. The world-line-based for- 
mulations of quantum theories are actively developed not 
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only due to the wide range of applications to the standard 
perturbative QFT, but also from the point of view of var- 
ious string theories (for recent review, see [20] and Refs. 
therein). 

The very important feature of this approach is that it 
does not refer directly to the standard perturbative tech- 
niques what allows one to avoid diagrammatic calculations 
which use to be very difficult and involved in non-Abclian 
gauge theories. Therefore, it is equally suitable for the 
perturbative as well as nonperturbative calculations [22]. 
The nonperturbative calculations within the the world- 
lincs approach could be compared with the results ob- 
tained in lattice QCD and other nonperturbative frame- 
works. The quark-vector boson vertex, being one of the 
fundamental objects in the theory of strong interactions, 
is under active investigation nowadays. The resummation 
of the gluon radiative corrections to the qqg-veitex has a 
direct relation to the problem of IR behavior, chiral sym- 
metry breaking, and confinement in QCD (see, Ref. [21] 
and Refs. therein). In high-energy scattering, it can con- 
tribute, e.g., to the amplitude with two-gluon exchange in 
the i-channel in the color singlet state. 

The RG equation and corresponding gauge-dependent 
cusp anomalous dimension which defines the IR properties 
of the non-singlet quark form factor were studied for the 
first time in Ref. [12] in the leading logarithmic approx- 
imation within the standard diagrammatic approach. In 
the present paper, we make an attempt to generalize the 
path-integral formulation to colored quark-vector boson 
vertices and found in one-loop order the cusp anomalous 
dimension of the averaged path-ordered Wilson exponen- 
tial corresponding to the color non-singlet quark form fac- 
tor. We point out an erroneous sign of one of the terms of 
the anomalous dimension calculated in [12], and find and 
check the correct expression. Thus, by virtue of exponen- 
tiation, one concludes that the colored quark form factor 
demonstrates the enhancement at large Q 2 , in contrast 
to the colorless case, Eq. (1). The origin of this enhance- 
ment is in close connection with the structure of the gauge 
group and has the same nature as the similar behavior in 
the decays of the charged electroweak bosons into fermions 
[19]. 



2 Soft gluon resummation within the 
world-line approach 

In this Section, we generalize the world-line method for 
evaluation of the QCD amplitudes with resummed radia- 
tive gluon corrections developed in [18] to the case of the 
quark-vector boson vertex. For this purpose, let us con- 
sider the color non-singlet form factor of a quark which 
can be extracted from the amplitude of the quasi-elastic 
(color changing) quark scattering (in the given kinemat- 
ics) in an external colored (gluon) gauge field: the quark 
on a mass shell comes from infinity, emits a hard gluon at 
the origin, changes the color, and goes away to infinity: 



F™ v [(pi- P2 ) 2 ;£] UifaWiWfa) 



(3) 



where pi <2 are the momenta of the in-coming and out- 
going quarks, and a, (i, j) are the gluon and quark color 
indices, respectively. This vertex is gauge non-invariant, 
therefore the dependence from the gauge parameter £ is 
included. We consider the covariant gauge with the fol- 
lowing gauge-fixing term in the Lagrangian: 



A 



i 

A 



(4) 



-the axial gauge C/^f 



while the other important case 

— X'/2(n- A) 2 — deserves a special attention. For the choice 
of kinematics with the large (longitudinal) momentum 
transfer between quarks (Sudakov regime), the interac- 
tions of quarks with external gluon, as well as the (external- 
)gluon-gluon and gluon-ghost interactions were shown to 
be IR-safe, thus they can be neglected in the considered 
case [11]. The kinematics is fixed by the small masses of 
the quarks and large squared transferred momentum: 



{piVi) — to 2 cosh x , (pip 2 ) > m 2 (5) 



or 



s = (pi +P2) 2 = 2m 2 (l 

-Q 2 = t = ( Pl -p 2 ) 2 = 

(s + 1) > (s - t) . 



-coshx) , 
2m 2 (l — coshx) , 



(6) 



According to the Feynman rules, the color structure of the 
vertex function (3) is determined by the matrix tfj in the 
fundamental representation of the Lie algebra of the color 
gauge group SU(N C ). 

Within the world-line formalism, the two-point (Eu- 
clidean) fermionic Green function can be written in terms 
of the Polyakov path-integral representation [17,18]: 



Gij(x,y) 



*i 3 

-i dr e~ m T / Vx{t') 

J J Xq—X 



m - -7 ■ x{t) 



7>cxp Q J dr' (TftvUft^j ■ cxp J dr' i 2 (r')^ 



•"Pcxp 



Here 



I dxpAp 
Jo 

~2 (*^ \x&v 



(x(r')) 



(7) 



is the Polyakov spin factor, 



&tiv — [linlAi are the quarks color indices, and r' 

is treated as the proper time of a fermion with mass m 
traveling from the point x at r' = to y at r' = r. Now 
we are interesting in the process where a fermion with 
initial momentum p\ and color i starts its way at point 
x, passes through a point z where it experiences both the 
large momentum transfer Q = p\ — p 2 , Q 2 < and the 
color changing described by the matrix T-j, and comes 
finally to point y with the momentum p 2 and color j. This 
process can be studied in terms of the three-point vertex 
function 



V*{x,y,z) = G ii ,(x,z)r$ j ,G j . j {z,y) , 



(8) 
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where the vertex matrix consists of a Dirac component 
r M = 7^,757^... and a color component T H = 5m, t^: 
rfa = <g) . For simplicity, we will consider the last 
case which corresponds to the quark-gluon vertex, while 
the color singlet quark-photon result is restored easily. 
Therefore, in the momentum space, the vertex function 
reads 



[PUP2) = 



--.(2) 



Hi 



■K-Pexp 



i5 / dx M 



,(2) 



(9) 



(1 2) 

where the sum over all possible trajectories C k / of two 
quarks is assumed. The functions [Ck,i] accumulate 
information about quark propagators. In this language, 
the extraction of the soft (long-distance) part can be per- 
formed by choosing a special set of paths with simple ge- 
ometry for ordered exponentials in Eq. (9). In our case, 
one should take for this purpose an angle with semi-infinite 
sides that represent the classical trajectories of the quarks 
[16, 18]. Then, the UV cutoff which arises in evaluation of 
the "soft" integrals is identified with the IR cutoff of the 
"hard" part. After rcnormalization of the soft exponen- 
tials (for details, see Refs. [15]), one can find the so-called 
cusp anomalous dimension which determines the large-Q 2 
asymptotic of the form factor 

F 1lV [q2] = 

(" ^2 ^ 7T + NLO ) - (10) 

By virtue of the Eq. (9), it is convenient to express the 
latter in terms of the vacuum average of the two path- 
ordered Wilson exponentials: 

T* F qqV [Q 2 ;£] = (o|T{w^.,W^.}|o} . (11) 

Taking into account that Eq. (11) contains, in general, 
both UV and IR divergences, we use scales /i 2 and A 2 as 
the UV and IR regulator, respectively. For the on-mass- 
shell quarks, their trajectories are the semi-infinitcly ex- 
tended paths, and can be parameterized as: 



In: 



re [-oo,0] , vl=pl/m, (12) 

Out: y v = v 2 a , a6 [0, +oo] , v 2 =p 2 /m. (13) 
Thus, the path-ordered exponentials can be written as 



Wu> = V exp 



and 



Wj/j = V exp 



ig t a vl f da A^a) 

J —GO 

ig fvl / da A%v 2 a) 
Jo 



(14) 



(15) 



The non-zero contributions (up to 0(g 2 ) terms) to F qqg 
(11) stem from the terms: 



W = 6 ii ,5 rj % j ,=1?, 



W 



(i) 

LO 



(16) 



pOO pOO 

/ da da' 6{a-a') D^, [v\a-a')] , (17) 
Jo Jo 

(2) 



W 



LO 



,2 

2 hj 



/0 rO 
da / da' 9(a - a') [v 2 (a - a')] , (18) 

-oo J — oo 

W LO = ~Y( CF ~ T") ■ 

pOO pOO 

/ dr da D^(v 1 T + v 2 a) . 
Jo Jo 



(19) 



Here the free gluon propagator is 

'o 



TA«(x)A?(y) 0) - Vf v {x - y) = S^D^x - y) , 

(20) 

and the following relations have been used: 
N 2 



*fc»*3 = —hr T — fat = C f 5m , 



2N r 



5u5jk - —5ij5ki 



(21) 



3 Leading order contributions for an arbitrary 
gauge field 

First, we evaluate the vacuum averaged Wilson integral 
Eq. (11) in n-dimensional space-time, for an arbitrary 
gauge field which can be of any origin, for instance — 
nonperturbativc. Then the general results will be applied 
in the case of usual perturbative gluon field. 

It is convenient to present the gauge field propagator 
D^ v (z) in the form [7]: 

D^(z) = g^dpOPD^z 2 ) - d^D 2 (z 2 ) . (22) 

Calculation of the integrals (17, 19) requires the following 
expressions: 

d p d p = 2nd + Az 2 d 2 , = 2g^d + Az^z v d 2 , d = d z 2 . 

(23) 

The scalar products are: 

= (<7 - °' f , v^vlg^ = coshx , 
v 1 ll v 2 z^ L z v = (vV + v 2 t) 2 coshx — ctt sinh 2 \ ■ (24) 

In order to control UV singularity, the transverse space- 
like separation b = (0ii,6j_) between two points at the 
integration path is introduced: 



3 3 



(a - a') 2 - (a - a') 2 - b 2 , 
{v 1 a + v 2 T) 2 ^{v 1 a + v 2 T) 2 -b 2 



(25) 



4 



Igor O. Cherednikov: Soft gluon resummation in quark-vector boson vertex at high energy 



The path-ordered integrals from Eqs. (17,19) can be 
evaluated using the following basic integrals: 



/ dada' e- a< -°-°"> 2 = -— , 
Jo 2a 



dada' aa' e'^-^ 



12a 2 



I 



dadr e -*(° 2 +T 2 +2vr cosh x ) = 



1 X 



2a sinh \ 

/ dadr ar e -«(- 2 +- 2 + 2 ^cosh x ) = 
Jo 

_ xcothx - 1 
4a 2 sinh 2 \ 



(26) 



Then, applying the Laplace transform to the invariant 
functions D^u), and it's derivatives over u — z 2 : 

/>oo 

D\ k} {u) = {-) k / da a k e' au Di(a) , (27) 
Jo 

one obtains the following formula: 

r°° i 

/ dada' D'{u) = -D(-b 2 ) , 
Jo 2 
f°° 1 

/ dada' D"{u) = -D'(-b 2 ) , (28) 
Jo 2 
f°° 1 

/ dada' aa' D" (u) = —D(-b 2 ) , (29) 
Jo 12 



Jo 



dada' (a - a') 2 D"(u) = --D(-b 2 ) (30) 



for u = (a — a') 2 - b 2 , and 

/>oo 

/ dadr D'(u) : 
Jo 



D{-b 2 ) 



POO 

/ dadr aTD"{u) 
Jo 



2 sinhx 
Xcothx - 1 



D(~b 2 ) , 



'0 

/•OQ 



4 sinh x 

/ dadr {v 1 a + v 2 T) 2 D"{u) 
Jo V > y > 2sinh X 

/•OO 



1 ^(-b 2 ) 



2 sinh x 



(31) 



for u = (uV + w 2 r) 2 - 6 2 . 

Thus, one finds the general expressions in n dimensions 
and arbitrary covariant gauge, with the gauge field two- 
point correlator expressed like the Eq. (22): 



W (l) = W (2) = _ t a 9_ Cf 



up to 0(g 4 ) order terms. Let us emphasize, that the equa- 
tions (32, 33) derived above are valid for any gauge field 
in the adjoint representation of the SU(N C ) color group. 
Therefore, these results can be used for evaluation of the 
nonpcrturbative, e.g., instanton, contributions to the vac- 
uum averaged Wilson integrals in the Eq. (11). 



4 One-loop perturbative contribution 

Now let us apply this general result in the particular case 
of the perturbative gluon field. Here and in what follows, 
the dimensional regularization is used with n = 4 — 2e, 
e < in order to regulate the IR-divergent terms in the 
integrals, respecting, in the same time, the gauge invari- 
ance. The gluon propagator in the coordinate space reads: 



r d n 

J (27T 



k e 



J [2w) n k 2 + iO 

r(i-£) 



(z 2 - iO) 



l-e 



k 2 +M 

r(-e) 
(z 2 - ioy £ 



(34) 



The perturbative one-loop invariant functions in Eq. (22) 
then read: 



Di(z 2 ) = 
D2(z 2 ) = 



- - ' r{1 ~ £) (-nX 2 z 2 ) e , 



167T 2 £ 

-iD.iz 2 ) . 



(35) 



Now the Eqs. (32, 33) contain the IR singularities at e — > 
0. Performing the standard rcnormalization procedure within 
the MS scheme, described in detail in Rcfs. [15, 14, 16], one 
finds 



and 



LO 
.a a $ s~i 



a. s ( „ C_a 
2 



i ^ 1 " 



(36) 



Xcothx 



In 



A 2 



(37) 



Here the UV- normalization point is taken /i 2 = 4fo~ 2 . 

Then one needs to combine the expression for the one- 
loop contribution to the form factor from Eqs. (16, 36, 
37): 



2 / 2 

[(n - 2)D 1 (-b 2 ) + 2b 2 D[(-b 2 ) + D 2 (-b 2 )] , (32) Flo&O = W ° + 2W lo ^ 



+ 



W^ 2 \x)=t%9 2 (C F -^ 

■ [ X cothx((n-2)D 1 {-b 2 ) + 
2b 2 D[(-b 2 )) +D 2 (-b 2 )] 



+W W(a s ,^,x)+0(a 2 s ). 



(38) 



(33) 



The high-energy asymptotic behavior of the form fac- 
tor is determined by the (gauge-dependent) cusp anoma- 
lous dimension which stems from the rcnormalization of 
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the Wilson integral (11) [11,12,16]: 



— 2 r^ sp [a s {^)- x ] , 

with the one-loop functions 

,da s 



In F N »(Q 2 ) = 



(39) 



0(a.) = ^ 



d^t 2 



-^-a 2 s + 0(a% ft, 



11 



-N r . 



(40) 



and (see Ref. [23]) 

4tt £ 



+ 



13 



+ 0{al) 



(41) 



Thus, the calculation of the cusp anomalous dimension 
(39) yields in one-loop approximation: 



- cusp 
7T 



C A \ C a 
Cf - — I xcothx + — £ 



0(^X42) 



where the Casimir operator of SU (N c ) in the adjoint rep- 
resentation is used: Ca = N c . Note, that in the Ref. [12], 
the term Ca£/4 in Eq. (42) had the sign minus. How- 
ever, our expression allows to reproduce the color singlet 
case by means of the replacement of the color factors: 
(Cf — Ca/2) — > Cf, thus obtainig the gauge invariant 
result (2): 



^ [C FX coth x + (C F -C f )£- C f ] + 0(a 2 s ) 

TT 

" 1 (xcothx - 1) + 0(a 2 s ) , 



2iV n 



(43) 



yielding the well known Sudakov suppression. The latter 
can be considered as a convenient test of the calculations 
and confirms the accuracy of our result, Eq. (42). 

In order to find the large-Q 2 asymptotic, we take into 
account the limit of the large scattering angle: 



Q 2 

Xcothx oc In — 



oo 



(44) 



and find that, in this limit, the anomalous dimension is 
linear in InQ 2 while the gauge-dependent term can be 
neglected: 



Note that for the color group SU(3), this anomalous di- 
mension (46) has the opposite sign compared to the color 
singlet case Eq. (2): this is a direct consequence of the 
algebra of gauge group generators. Therefore, the lead- 
ing (double-logarithmic) behavior of the non-singlet form 
factor is given by 

F 119 [Q2] ipqqg [Q2] = 



exp 



exp 



- C F 



rQ dfi ^Q 2 ^ a s ((i) 
, Q 2 In n TT 



NLO 



PoN c 



In 



Q 2 ^ \nQ 2 /A 2 QCD 

^ ^Qo/ A QCD 



A 2 



+ NLO 



(47) 



that is the increasing function of Q 2 . The dependence from 
the gauge-fixing parameter £ drops out of the leading log- 
arithmic expression for r™ p , Eq. (46), and yields no in- 
fluence on the main asymptotics. 



5 Discussion and Conclusions 

We derived the general formula for the vacuum averaged 
Wilson integral (11) in the g 2 accuracy. This result Eqs. 
(32, 33) is found in n dimensions, in an arbitrary covariant 
gauge, and is valid as well for any SU (N c ) gauge field. In 
the case of a nonperturbative field A^ p <~ g -1 , this result 
corresponds to the leading (so-called "weak field") order 
of expansion in field strength, while the dependence on 
the coupling g drops out. 

By using the Wilson integrals techniques, it has been 
found that the cusp anomalous dimension for the non- 
singlet quark form factor Eq. (46) has a negative sign — in 
contrast to the singlet case — what leads, in the large-Q 2 
regime, to the enhancement, rather than suppression, of 
the contribution due to the resumed soft gluon radiative 
corrections, Eq. (47). The dependence from the gauge- 
fixing parameter £ is shown to be of the order of In Q 2 — 
not surviving in the asymptotical regime. 

It is necessary to emphasize that we work here in the 
covariant gauge, whereas the case of the axial gauge is 
quite important and requires additional analysis (see, e.g., 
Rcfs. [24]). Another important version of the problem 
the case of the q/f-shell quarks with different masses (ap- 
plicable, e.g., to flavor changing processes) is more tech- 
nically involved and will be reported in the forthcoming 
work. 

Another point which should be noted is that the quark- 
gluon vertex is a colored gauge-dependent object and could 
not be an observable quantity. In computations of the re- 
alistic processes when this vertex is inserted into the dia- 
grams at the partonic level, the role of this enhancement 
can be reduced due to the emission of gluons, as well as 
due to the convolutions of the color indices. 



HSSp [a.; X ] = r™% [a s ] In ^ + O(ln Q 2 ) , (45) Acknowledgments 



-^cusp [ a s] — 

TT 



C> ^1 <)• (46) 



The author thanks A.E. Dorokhov, L. Magnea and N.G. 
Stefanis for careful reading of the manuscript, fruitful dis- 



6 



Igor O. Cherednikov: Soft gluon resummation in quark-vector boson vertex at high energy 



cussions and valuable remarks. The conversations with 
E.A. Kuraev and O.V. Teryaev on a number of topics con- 
cerning this work were very useful and inspiring. The work 
is partially supported by RFBR (Grant Nos. 04-02-16445, 
03-02-17291), Russian Federation President's Grant 1450- 
2003-2. The hospitality and financial support of the INFN 
Section in Torino (where this work was started), and the 
High Energy Section at ICTP in Trieste (where it was 
completed) is gratefully acknowledged. 



References 

1. S. Catani et al, hep - ph/0005025; S. Catani, 
hep - ph/0005233. 

2. P. Ciafaloni, D. Comelli, Phys. Lett. B476 (2000) 49; V.S. 
Fadin, L.N. Lipatov, A.D. Martin, M. Melles, Phys. Rev. 
D61 (2000) 094002; J.H. Kuhn, A. A. Pcnin, V.A. Smirnov, 
Eur. Phys. J. C17 (2000) 97. 

3. Yu. Dokshitzer, D. Dyakonov, S. Troyan, Phys. Reports 58 

(1980) 269; A. H. Mueller, Phys. Reports 73 (1981) 237. 

4. G. Parisi, Phys. Lett. B90 (1980) 295; G. Curci, M. Greco, 
Phys. Lett. B92 (1980) 175; H.n. Li, G. Sterman, Nucl. 
Phys. B381 (1992) 129; L. Magnea, G. Sterman, Phys. 
Rev. D42 (1990) 4222; L. Magnea, Nucl. Phys. B593 
(2001) 269; 

5. R. Pctronzio, S. Simula, G. Ricco, Phys. Rev. D67 (2003) 
094004; S. Simula, Phys. Lett. B574 (2003) 189; M. Os- 
ipenko et al. [CLAS Collaboration], Phys. Rev. D67 (2003) 
092001; M. Osipenko et al, hep - ex/0309052. 

6. N. Kochelev, Phys. Lett. B565 (2003) 131. 

7. A. Dorokhov, I. Cherednikov, Phys. Rev. D66 (2002) 
074009; Phys. Rev. D67 (2003) 114017; Annals Phys. 
(NY) 314 (2004) 321; I. Cherednikov, Surv. High Energy 
Phys. 18 (2003) 205. 

8. V. V. Sudakov, Sov. Phys. JETP 3 (1956) 65; [Zh. Eksp. 
Teor. Fiz. 30 (1956) 87]. 

9. E. A. Kuraev, V. S. Fadin, Yad. Fiz. 27 (1978) 1107; B.I. 
Ermolaev, V.S. Fadin, JETP Lett. 33 (1981) 269; [Pisma 
Zh. Eksp. Teor. Fiz. 33 (1981) 285]; B.I. Ermolaev, S.I. 
Troyan, Nucl. Phys. B590 (2000) 521. 

10. R. Jackiw, Annals Phys. (NY) 48 (1968) 292; P. M. Fish- 
bane, J. D. Sullivan, Phys. Rev. D4 (1971) 458; J. M. 
Cornwall, G. Tiktopoulos, Phys. Rev. D13 (1976) 3370; 
J. Frenkel, J. C. Taylor, Nucl. Phys. B116 (1976) 185; E. 
C. Poggio, G. Pollack, Phys. Lett. B71 (1977) 135; A. H. 
Mueller, Phys. Rev. D20 (1979) 2037; J. C. Collins, Phys. 
Rev. D22 (1980) 1478; A. Sen, Phys. Rev. D24 (1981) 
3281. 

11. J.J. Carazzone, E.C. Poggio, H.R. Quinn, Phys. Rev. Dll 
(1975) 2286; [Erratum-ibid. D12 (1975) 3368]. 

12. Yu. A. Sitenko, Sov. J. Nucl. Phys. 27 (1978) 583; [Yad. 
Fiz. 27 (1978) 1098 (in Russian)]. 

13. A. Bassetto, M. Ciafaloni, G. Marchesini, Phys. Reports 
100 (1983) 201; O. Nachtmann, Annals Phys. (N.Y.) 209 
(1991) 436. 

14. Yu.M. Makeenko, A. A. Migdal, Phys. Lett. B88 (1979) 
135; Nucl. Phys. B188 (1981) 269; A. Polyakov, Phys. 
Lett. B82 (1979) 247; Nucl. Phys. B164 (1980) 171. 

15. V. S. Dotsenko, S. N. Vergeles, Nucl. Phys. B169 (1980) 
527; R. A. Brandt, F. Neri, M.-A. Sato, Phys. Rev. D24 

(1981) 879; R. A. Brandt, A. Gocksch, M. A. Sato, F. Neri, 



Phys. Rev. D26 (1982) 3611. J.G.M. Gatheral, Phys. Lett. 
B133 (1983) 90; J. Frenkel, J.C. Taylor, Nucl. Phys. B246 
(1984) 231. 

16. G. Korchemsky, Phys. Lett. B217 (1989) 330; G. Ko- 
rchemsky, A. Radyushkin, Sov. J. Nucl. Phys. 45 (1987) 
910; G. Korchemsky, A. Radyushkin, Nucl. Phys. B283 
(1987) 342; G. Korchemsky, Phys. Lett. B220 (1989) 629. 

17. A.I. Karanikas, C.N. Ktorides, Phys. Rev. D52 (1995) 
5883; A.I. Karanikas, C.N. Ktorides, N.G Stefanis, Phys. 
Rev. D52 (1995) 5898. 

18. GC. Gellas, A.I. Karanikas, C.N. Ktorides, N.G. Ste- 
fanis, Phys. Lett. B412 (1997) 95; N.G. Stefanis, 
hep - ph/9811262; A.I. Karanikas, C.N. Ktorides, N.G. 
Stefanis, Eur. Phys. J. C26 (2003) 445. 

19. B.I. Ermolaev, hep - ph/0211302; B.I. Ermolaev, S.M. 
Oliveira, S.I. Troyan, Phys. Rev. D66 (2002) 114018; 
A.Barroso, B.I. Ermolaev, hep - th/0112216. 

20. C. Schubert, Phys. Reptorts 355 (2001) 73. 

21. S.J. Stainsby, R.T. Cahill, Mod. Phys. Lett. A9 (1994) 
3551; A.I. Davydychev, P. Osland, L. Saks, Phys. Rev. 
D63 (2001) 014022; A. Bender, W. Detmold, CD. 
Roberts, A.W. Thomas, Phys. Rev. C65 (2002) 065203; 
J.I. Skullerud, P.O. Bowman, A. Kizilersu, D.B. Leinwe- 
ber, A.G Williams, JHEP 0304, 047 (2003). 

22. Yu.A. Simonov, J.A. Tjon, Annals Phys. 300 (2002) 54; 
A.V. Belitsky, A.S. Gorsky, GP. Korchemsky, Nucl. Phys. 
B667 (2003) 3. 

23. O.V. Tarasov, A. A. Vladimirov, Sov. J. Nucl. Phys. 25 
(1977) 585; [Yad. Fiz. 25 (1977) 1104]. 

24. G.T. Bodwin, S.J. Brodsky, GP. Lepage, Phys. Rev. Lett. 
47 (1981) 1799; R. Doria, J. Frenkel, J.C. Taylor, Nucl. 
Phys. B168 (1980) 93. 



